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We apply the extended (statistically-consistent, SGA) Gutzwiller-type approach to the periodic 
Anderson model (PAM) in an applied magnetic field and in the strong correlation limit. The finite- 
U corrections are included systematically by transforming PAM into the form with Kondo-type 
interaction and residual hybridization, appearing both at the same time. This effective Hamiltonian 
represents the essence of Anderson-Kondo lattice model. We show that in ferromagnetic phases the 
low-energy single-particle states are strongly affected by the presence of the applied magnetic field. 
We also find that for large values of hybridization strength the system enters the so-called locked 
heavy fermion state. In this state the chemical potential lies in the majority-spin hybridization gap 
and as a consequence, the system evolution is insensitive to further increase of the applied field. 
However, for a sufficiently strong magnetic field, the system transforms from the locked state to 
the fully spin-polarized phase. This is accompanied by a metamagnetic transition, as well as by 
drastic reduction of the effective mass of quasiparticles. In particular, we observe a reduction of 
effective mass enhancement in the majority-spin subband by as much as 20% in the fully polarized 
state. The findings are consistent with experimental results for CexLai-ccBe compounds. The mass 
enhancement for the spin-minority electrons may also diminish with the increasing field, unlike for 
the quasiparticles states in a single narrow band in the same limit of strong correlations. 

PACS numbers: 71.27. +a, 71.10.Fd, 71.28.+d, 75.30.Mb 



I. INTRODUCTION 

Heavy fermion (HF) metals are customarily regarded 
as a quantum liquid composed of strongly correlated 
(quasi-atomic) and heavy f-electron quasiparticle states 
resulting from hybridization with conduction (band) elec- 
trons. The starting f-electrons are associated with the 
presence of rare-earth or actinide ions (ions with valency 
close to integer valence, e.g. Ce +3_5 , with 5 <C 1). 
In other words, the f-f on-site Coulomb repulsion is 
strong and the f-electrons behave often in a solid as 
localized (or almost localized) magnetic moments even 
in the presence of the interband hybridization. More- 
over, those heavy hybrid f-electrons are coupled to the 
conduction (d-s) electrons via antiferromagnetic Kondo 
inter act ioni. Also, many of the cerium or uranium 
HF compounds order antiferromagnetically. Nonethe- 
less, there is also a constantly growing group of heavy 
fermion compounds which show a ferromagnetic order- 
ing, e.g. URu 2 - a ,Re :c Si2 for 0.3 < x < l.OS, UIr 2 Zn 2a 2, 
CeSii.8i 4 ,YbRh 2 Si2^. Therefore, it is important to un- 
derstand the complexity of their phase diagram within 
a single theoretical framework such as Anderson— or 
Kondo-lattice models 7 . In this respect, we use here a 
mixed Anderson-Kondo lattice picture derived from PAM 
in a systematic manner. 

One of other characteristic properties of the heavy 
fermion compounds is the metamagnetic behavior or a 
sudden magnetization increase (metamagnetism) at a 
critical value of the applied magnetic field, observed in 
HF compounds such as CeRu 2 Si2^, CeCoGe^, UPd2Al3, 
URu 2 Si2, and UPt^. This feature of HF electron sys- 



tems makes the applied field a useful tool of probing 
the strong-correlation behavior. For example, a suffi- 
ciently strong applied field may lead to a destruction of 
the HF state of those strongly correlated fermions and, 
as a consequence, force a localization (freezing) of the 
heavy, predominantly f-electron, quantum liquid. This 
type of transition signalled by a metamagnetic transition 
(MMT), was the subject of extensive studies during the 
past decade, but still, because of the rich variety and 
complicated structure of their phase diagrams, the uni- 
versal nature of the MMT is hard to establish. One of 
the puzzles is that the metamagnetic behavior can take 
place starting from the ferromagnetic state. 

In this paper we present a scenario in which the heavy 
itinerant f-electron quasiparticles localize at MMT, with 
a concomitant suppression of the HF state at that critical 
field. To account for this interplay between the localized- 
and the itinerant-type behavior we start from the pe- 
riodic Anderson model (PAM) in an external magnetic 
field included via the Zeeman term. We also use the 
large, but finite U version of PAM formulated befor e - 11 ! 12 , 
in which both residual hybridization and the Kondo- 
type coupling are included in a systematic manner. We 
have used this model befor o 13 i 14 , which can be called 
the Anderson-Kondo lattice model, to discuss the so- 
called hybrid (Kondo-type) pairing in heavy-fermion sys- 
tems. Additionally we utilize the so-called statistically- 
consitent extension of the Gutzwiller-type approac h 15 ! 16 . 
Such formulation allows for a consistent mean-field dis- 
cussion of correlated states, here in the applied magnetic 
field. With this paper we are starting a series with a 
detailed analysis of various magnetic, as well as of the 
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magnetic superconducting states. 

The paper is organized as follows. In Sec. El we sum- 
marize our theoretical model. In Sec. IIIII we discuss 
the so-called statistically consistent mean-field approach 
(SGA) and provide numerical results, as well as apply 
them to the analysis of Ce^Lai-^Be properties. We sum- 
marize our results and provide a brief outlook in Sec. HVl 
The Appendices IX1IO provide some details of our analytic 
part of the approach. 



II. MODEL AND STATISTIC ALLY-CONSITENT 
MEAN FIELD THEORY 

A. Anderson-Kondo lattice Hamiltonian 

In modeling the fluctuating- valence and heavy fermion 
systems one usually starts from the periodic Anderson 
model (PAM) 17 , which in the real space (Wannier) rep- 
resentation reads 
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FIG. 1: (Color online) Schematic representation of the first 

(a) and second-order (b-d) hopping processes. The hopping 
label 2 and 2' in (c) are alternative processes. The process 

(b) and (c) lead to real-space pairing, whereas the three-step 
process (d) leads to itineracy of originally atomic f electrons 
with effective hopping tij . All the processes contribute to the 
dynamics of heavy quasiparticles with renormalized charac- 
teristics. The virtual process (b) leads to the Kondo interac- 
tion, whereas the process (c) introduces hybrid-pair hopping. 
The effective f-f hopping (d) results from residual hybridiza- 
tion in the narrow-band limit (for detailed discussion of this 
last picture see the next Subsection III B|) . 



with Ni a = fj a fia being the number of originally atomic 
(f) electrons and h ma = c^Cma the number of conduct- 
ing (c) electrons. The consecutive terms are: the first 
represents the hoping (band) energy of c electrons, the 
second - bare f atomic level position (with respect to the 
atomic level of c electrons, ey = ey — e c ), the third - inter- 
atomic Coulomb interaction among the f-electrons (the 
Hubbard term), the fourth the f-c hybridization (with 
the amplitude V^ m , taken in the subsequent analysis in 
the intraatomic form, Vi m = Si m V), and the next term 
represents the subtraction of the chemical potential part, 
as we will work in the grand canonical scheme. The last 
two terms are the Zeeman energy in the applied magnetic 
field H of f and c electrons, respectively. 

In the context of heavy fermion system this Hamil- 
tonian has usually been taken in the U = oo limit, as 
this parameter represents the highest energy scale in 
the system. Typically for Ce systems: U = 5 — 6eV, 
e f = -1 -2eV, V < -0.5eV, W = 2z\t\ ~ 2eV-±2. 
While |F|/?7 <C 1, U is definitely finite and therefore, 
the finite-U corrections should be taken into account, as 
e.g. the residual Kondo exchange interaction of the mag- 
nitude J K = 2V 2 /(e f + U) - O.W arises 12 and will 
influence the nature of the magnetic ground state in a 
decisive manner. However, to account for the Kondo in- 
teraction systematically, as well as to allow for itineracy 
of f-electrons at the same time, a direct application of ei- 
ther the Schrieffer-Wolf£i2 transformation or by starting 



from the periodic Kondo model 7 (i.e. neglecting explic- 
itly the hybridization term apart from the presence in the 
expression ~ J^), may not be realistic. At least, one has 
to try to reach such Kondo-lattice limit in a systematic 
manner, by e.g. including explicitly both the nonzero 
renormalized hybridization and the nonvanishing Kondo 
coupling in the large-U limit at the same time. 

For this purpose, in direct analogy to the correspond- 
ing projection set originally for the Hubbard model 2 ^, we 
divide the f-c hybridization into the following two parts 
according to the prescription 

= {l-N is )flc ma + N l9 fic ma . (2) 

The first term represents the projected hybridization of 
the states \fi a ) and |c mcr ) under the proviso that the f- 
state \fia) is not occupied, whereas the second represents 
the part of the f-c quantum-mechanical mixing when the 
admixture of a doubly-occupied state with the opposite 
spins, \fi a ) and \fi a ), is admissible. The various hopping 
processes involving the original (bare) atomic f electrons 
and the resulting occupation-resolved hybridization pro- 
cesses are depicted in Fig. Q] (for their detailed meaning 
see below). Since the heavy fermion physics is related to 
the limit: \ef\ > V and \ef\ <C e / + U (with e f = e f — ji) , 
the first term in (|2]) corresponds to the low-energy mixing 
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processes (represent a real f-c mixing, Fig. D^b)), whereas 
the processes represented by the second term are realized 
only via virtual (high-energy processes, Fig. QJc,d)) and 
are accounted for in the second-order in V/(U + e/) (the 
first nontrivial order). The latter processes are removed 
from the original Hamiltonian (pp) via canonical transfor- 
mation proposed earlie r 12 ! 21 . The original Hamiltonian 
is then transformed to the effective Hamiltonian of the 
form 
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where t rnn = t mn — //, cj = ej — jjl. Also, in this Hamilto- 
nian we have projected out completely the double occu- 
pancies of f states, as well as have defined the following 
projected quantities for f electrons Oif a = Ni a (l — iV^), 

% = £ CT Oifa, s, = (Sf, Sf) = [flu, l/2(% t - P ifi )], 
and the corresponding (non-projected) quantities for the 
conducting band electrons are n mcr = cj^c^, n m = 
Ea ™ma and s m (s^, sf) = [c^Cma, l/2(ra mt - n m ;)]. 
Additionally, as we assume that g c = = g = 2, 
the applied field is defined as ft, = ^g/isH. Finally, 



we have added an effective 4-th order term 



^2 m (VimVjm) 4: / (ej+L 7 ") 3 , which will not be written up ex- 
plicitly here, as we concentrate only on the second-order 
effects. The Gutzwiller projector V = JJ^l - N^N^) 
eliminates explicitly the double f-occupancies (it can be 
dropped as the operators inside V{...)V are already in the 
projected form). The first term in (|3]) contains c-c hoping 
(see Fig. Ufa)) renormalized by the virtual- hopping pro- 
cess, the second term is effective hopping with the spin- 
flip exchange between c and f electrons (see Fig. Gib)), 
the third part and the fourth - the f-level energy and 
projected hybridization, while the fifth - the effective f-c 
antiferromagnetic Kondo coupling. The minimal model 
coming out of that procedure is to neglect an effective f-f 
exchange and the renormalization of the c — c hopping. 
In result, we start from a simplified effective Hamiltonian 
which we rewrite in the explicit form 
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We see that we have simultaneously presented a resid- 
ual (projected) hybridization term and the Kondo kinetic 
exchange interaction, with the exchange integral Jk = 
Ji m = 2\Vi m \ 2 /(U + which differs from that obtained 
with the help of the Schrieffer- Wolff transformation (as it 
contains only high-energy virtual- hopping processes, cf. 
Fig. DIb)). The appearance of a number of magnetic 
phases arises from a competition/cooperation of the pro- 
jected hybridization and antiferromagnetic Kondo inter- 
action. 



B. A general comment: Anderson, Kondo or mixed 
Anderson-Kondo lattice representation? 

It may be useful at this point to characterize briefly 
the subtle difference between the present formulation 
and the original Anderson- or Kondo-lattice models. 
First, our canonical transformation differs from the 
standard Schrieffer- Wolff transformation. Namely, in 
the Schrieffer- Wolffl^ (and related Coqblin-Schrieffer— ) 
transformations the whole hybridization term is trans- 
formed out and replaced by the Kondo interaction. Here, 
only the second part of ([2]) represents virtual processes 
(cf. Fig. DIb) ), whereas the first term of (|2j), cannot be 
transformed out, since V ~ ef and hence the correspond- 
ing processes lead to the itineracy of f-electrons (cf. Fig. 
Hfd) for an illustration of an effective f-f hopping). This 
means that the f-electrons are not strictly localized, as in 
most of the HF systems they are not. As a consequence, 
the total number of f particles is not conserved, only the 
total number N e = Y^ia ^ia + Y^ma ™rna of all parti- 
cles is a conserved quantity. The last statement amounts 
to saying that the physics of Hamiltonian (j3j) and (4) is 
contained in that coming of (pQ), but in general, is not 
equivalent to that of true Kondo-lattice model, i.e. when 
the third term in (4) ~ V is absent. This is similar to 
the difference between the physics coming from the An- 
derson impurity model as compared to that coming from 
the Kondo model; the later representing the asymptotic 
limit of the former—. 

In principle, one could say that all the relevant physics 
is fully contained in the general form (pQ). Why then 
introducing (J3j) or (4) which represent its particular lim- 
iting forms? This is the first basic question. The answer 



to this question is as follows. The form (pQ) represents 
indeed a general hybridized two-orbital system with the 
short-range (Hubbard) interaction (apart from the ab- 
sence of Falicov-Kimball term ~ X^<im> UimNifim, the 
role of the missing Falicov-Kimball term is discussed in 
Appendix [A]) . No exact solution is available in the lat- 
tice situation of dimensionality D » 1. Therefore, we 
have to resort to approximations and it is usually prof- 
itable, and even physically plausible, to take into account 
the principal interorbital/intersite exchange interaction 
explicitly first and carry out a simplified (e.g. mean- 
field like) analysis on such an effective Hamiltonian sub- 
sequently, with the most relevant correlations included, 
even in the approximate manner in the mean-field type 
approximation. What is equally important, the form of 
the effective Hamiltonian, by containing physical (low- 
energy) interactions, helps in selecting relevant order pa- 
rameters and related to them mean fields in nontrivial 
cases, as discussed in detail below in one such particu- 
lar situation. Parenthetically, the analysis is of the same 
type as that when transforming the Hubbard model into 
an effective t — J model and determining nontrivial many 
body ground states for the latter. Such Hamiltonian is 
applicable to the analysis of both normal and supercon- 
ducting states, the latter with real space pairing included. 

One has to note that for the Kondo-lattice Hamilto- 
nian, which has general form 
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in which the numbers of f-particles and c-particles are 
conserved separately, i.e. both Nf = ^2 io . Ni a and N c 



commute with %kl separately. Therefore, in 



discussing the HF or superconducting states based on this 
type of Hamiltonian^ one has to assume that those two 
global conservation laws are only approximately obeyed. 
It is better to use our effective Hamiltonian (|3j) or (4), 
on the expense of the simplicity though. 

A second basic question still remains and concerns the 
reduction of Hamiltonian (4) or (j3]) to (|5j) in a consistent 
manner. This is usually done by removing the residual 
hybridization term and completing the Schrieffer- Wolff 
transformation. However, in our mean-field analysis the 
true Kondo-lattice limit expressed by (|5j), is effectively 
achieved as a limit rif — >• 1 . This limit is indeed achieved 
when V/e — >• 0, as we show explicitly in the next Section 
in concrete situation. 



C. Efective mean- field Hamiltonian and 
statistically-consistent mean-field approximation 

In order to construct a mean field Hamiltonian we as- 
sume, that a proper ground variation state can be repre- 
sented in the form = < P|\E r o) 5 where indicates the 
correlated state and |^o) is an eigenstate of the effective 



single-particle Hamiltonian. Then the expectation value 
of any operator can be in principle calculated as 



(6) 
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In most cases, the above calculations represent a non 
trivial task, mostly because the projection operator is 
non-local and by applying the Wick theorem many 
terms appear in the process of Hartree-Fock type de- 
coupling. To carry out these evaluations an approxi- 
mation is necessary at this point. Here, we use the 
scheme proposed recently by Fukushima^ 5 - in the local- 
constraint version, which assumes that the average num- 
ber of correlated f-electrons at each site and with an 
arbitrary spin orientation is unchanged by the projec- 
tion, i.e. (Nia) — (A^icr)o = n ifa- I n the subsequent 
analysis, we use more general Gutzwiller-type projector, 

^ = Ui \ i t' t/2 \ i I u/2 ( 1 - Nilftit), which eliminates the 
double occupancies of f-electrons in real space and con- 
serves average number of f-electrons at each site before 
and after projection; the latter is accomplished by fix- 
ing the fugacity factors {A^}. In order to derive the 
explicit form of the fugacity factors, one has to calcu- 
late the density of a spin f-electron at site i, namely, 



. This can be done by 



(n z ^ 2 )o 

taking only the onsite contraction, as the intersite ones 
are much smaller. Then, / ^ i terms cancel out (in the 
numerator with these of the denominator). Therefore, 
the condition to determine the fugacity factors is given 
by (Nia) = which results in the ex- 

plicit expression for the fugacity factors, A^ = 25 • 

It should be noted that the uncorrelated c-electrons are 
not effected by the projection procedure, as it projects 
out only for the number of f-electrons. Though the 
original local-constraint scheme for calculating the cor- 
related averages was proposed for a single band Hubbard 
model 25 and was successfully applied in t — J model 2 ^, it 
can be easily extended to the periodic Anderson model in 
a straightforward manner (see Appendix iBl for details). 

In effect, this kind of Gutzwiller approximation yields 
the mean value of the Hamiltonian (4) 



A 



8i£ 



V 2 



V 2 



U + e f 



U + e f 



Si 



V 2 



U + e f ^ ^(1 -"/*)(! "»/*) 
ejrif — hfTrif — h c m c , 



1°% 



(7) 



5 



where A is the number of lattice sites and we have ad- 
ditionally assumed, that all the mean fields appearing in 
our model represent spatially homogeneous quantities, 
namely, i. e., 



nif = ^2 ^ ia ^ Umc = ^2 (^rna) n c , 

mi f = \ = m /' 

m mc = ^ ^ (crn m(7 ) ra c , 
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It should be noted that the hybridization term in ([7]) 
has the same renormalization factor as can be obtain in 
standard Gutzwiller approximatio n 27 ! 28 in the U = oc 
limit. 

Also, in order to formulate a completely correct mean- 
field approximation, we supply the mean-field Hamil- 
tonian ([7j) with constrains via the Lagrange-multiplier 
method. By doing so we ensure the statistical consis- 
tency of the present variational metho d 15 ! 26 of evaluat- 
ing the averages in the sense, that they also coincide 
with those calculated via a fully self-consistent proce- 
dure. This method provides results which are equiva- 
lent to those obtained in the saddle-point approximation 
within the slave-boson approach 2 ^. In result, the effective 
MF Hamiltonian takes the form 
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(9) 



where < m, n > means the summation over nearest 
neighboring pairs. Note that the field A/ m may play 
role of a molecular field induced by electronic correla- 
tions whereas Xf and A represent shifts of the chemical 
potential (the role of A proves to be minor). 

Now, Eq. ([9]) has a particularly simple k-space repre- 
sentation: 
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where n e = rtf + n c is the total number of electrons 
per site and the constant C includes all the wavevector 
independent quantities: C = (H) + A(16£?7 + A cm m c + 
^Scr r ^7cr + ^f n f + ^fm^f + An). For example, the 
band energy of c-electrons in a two-dimensional case is 
represented explicitly by e£ CT = — Arj (cos k x + cos k y ) — 
H — A — A cm <r and the energy of the localized f-electrons 



equals e } 



-fi — A — Xf — Xf m (j. It is easy to see that 



the eigenvalues of the Hamiltonian ([TO]) are then 
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where a = ±1 labels the quasiparticle band (antibonding 
and bonding components, respectively). The quasiparti- 
cles form two bands with renormalized characteristics. 

The mean-field values can be determined by min- 
imizing the free energy potential T = To — 
Xj3 ^kcra m (l + e~P E £°), i.e. by solving the following 
system of equation representing the necessary conditions 
for the minimum: 
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where equations with A^ = {A, Aj, Aj m , A cm , 77, r^} can 
be easily solved analytically, so that we are left with six 
equations which we solve numerically with the help of 
GSL (Gnu Scientific Library), on the lattice of A = 256 x 
256 sites, using the periodic boundary conditions, with 
the following values of the parameters: e/ = — 0.75W, 
t = 1/8 W, U = 5VK, and for sufficiently low temperature 
ksT/W = 5 • 10 -5 , where W is the bare conduction- 
electron bandwidth, that sets our energy scale (W = 1). 
Exemplary values obtained via the self-consistent and 
variational method for a ferromagnetic state and h = 
are displayed in Table HI The numerical accuracy is at 
least on the last digit of the quantities displayed in that 
Table. Also the numerical results are illustrated further 
on the example of a two-dimensional (square) lattice. 
Parenthetically, a detailed shape of the density of states 
in bare band is of secondary importance. 



III. RESULTS AND DISCUSSION 

A. Phase diagram and principal characteristics of 
selected phases 

A variety of magnetic phases appears in this model 
depending on the value of microscopic parameters and 
the type of long-range order obtained. We analyze in de- 
tail only the phases specified in Fig. [2] (bottom part) as 
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FIG. 2: (Color online) Phase diagram (top) on n e — \ V\ plane 
displaying various arrangements of quasiparticle spin sub- 
bands (bottom) relative to the chemical potential. In phase II 
(weakly ferromagnetic metallic - WFM) the lower majority- 
spin subband is fully occupied, whereas the minority-spin sub- 
band remains partially filled. The onset of region II is marked 
by the disappearance of the spin-up Fermi surface sheet, as 
the chemical potential is placed in the hybridization gap and 
the lower spin- up subband is completely filled. For the small 
values of hybridization, the majority-spin electrons fill the 
upper hybridized band, what results in a full polarization of 
the system (phase III - strongly feromagnetic metallic SMF). 
The case (IV) represents the paramagnetic Kondo insulator 
(PKI) and can be observed only for n e = 2.0. The phase 
diagram is drawn for U = 5, e/ = —0.75, and t — —1/8, all 
in units of the conduction electron bandwidth W. The differ- 
ence between various ferromagnetic phases is connected with 
the difference in the shape and the filling of various parts of 
the density-of- states, as shown at the bottom. 



TABLE I: Values of the chemical potential and the MF pa- 
rameters. The physical free energy F — T -V \in & for n e — 1.9 
and 1.75 is -1.04971 and -1.04175, respectively. We see that 
in this limit the f-electron occupation is close to unity and 
their magnetic moment (rrif) is decisively larger than that of 
the compensating Kondo cloud (ra c ). 
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well as provide the phase diagram on the n e — V plane 
(top part). For the case when the hybridization exists 
only between the electrons and holes at the same site, 
i.e. Vim — VSi m and is k- independent, we can precisely 
evaluate the effective mass of quasiparticles using analyt- 
ical results (see Appendix [Cl for details). It can be seen 
from Fig. [2] that, depending on the number of total elec- 
trons n e and the value of hybridization parameter |V|, 
we can achieve different states. In phase II, which cor- 
responds to the weakly ferromagnetic metal (WFM), the 
lower majority-spin subband is fully occupied, whereas 
the minority-spin subband remains partially filled. The 
onset of region II is marked by the disappearance of the 
spin- up Fermi-surface sheet, since the chemical potential 
is placed then in the hybridization gap and the lower spin- 
up subband is completely filled. If the lower majority- 
spin subband in phase II is fully occupied (contains one 
electron per site) then the minority-spin subband remains 
partially filled and contains n e — l electrons per site, what 
results in total magnetization m = (2 — n e )/2 and ap- 
proaches zero for n e — > 2 (cf. Fig. [3fa,c) and phase II 
in the Fig. [2]). On the other hand, for a small value of 
hybridization the majority-spin electrons in the phase III 
start filling the higher hybridized band, what results in 
an almost full polarization of the system (see Fig. [3fa,c)). 
As the magnetizations of f- and c-electrons, rrif and ra c , 
have opposite signs due to antiferromagnetic Kondo cou- 
pling between them, the maximum value of the total 
magnetization never exceeds 0.5/i# i.e. the maximum 
value in the localized limit. 

In the case of n e < 1.0 both f- and c-electrons have 
the same sign of magnetization, as a consequence of the 
fact that the effective hybridization, represented by the 
mean- field value of 7^ = (//^c^), is close to zero, as can 
be seen in Fig|4]b. In the case when n e > 1.0 and the bare 
hybridization \V\ < 0.3, the f-electrons start localizing, 
as can be seen by a rapid increase of rrif (see Fig. |3fa)) 
and the suppression of the hybridization between local- 
ized and itinerant electrons (in Fig. [4] see a recess of both 
7^ and 7^ in the region corresponding to the phase III 
(SFM)). The Kondo singlet correlation are destroyed to 
yield a strong ferromagnetism. Thus large magnetic mo- 
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FIG. 3: (Color online) The magnetic moments of f-electrons 
(a) and of conduction electrons (b) as function of bare hy- 
bridization magnitude \V\ and total number of electrons per 
site: n e = n c + n/, and (c) the number f-electron occupancy. 
Other parameters are the same as those listed previously. 



ments reflect the tendency to destroy the partial Kondo 
effect. Similar conclusions can be drawn from the S = 1 
underscreened Anderson lattice mode l 3Q i 31 . Here, the 
same effect is observed in orbit ally nongenerate PAM. 



B. Ferromagnetic phase: suppression of heavy 
quasiparticle masses 

The detailed phase diagram should encompass also an- 
tiferromagnetic phases 29 . In the present paper we shall 




FIG. 4: (Color online) The mean magnitude of hybridization 
correlation functions: 7 t = (f^c^) (a) and 7^ = (//^c^) 
(b) both, as a function of total number of electrons and bare 
hybridization parameter \V\. Other parameters are the same 
as those used previously. 



limit our discussion to the range of parameters, for which 
the ferromagnetic phases II and III are stable. The phase 
transition from the state with relatively small magnetic 
moments (phase II) to the phase with almost fully polar- 
ized spin-up subband (III) can be induced by reducing 
the hybridization strength |V|, as can be seen in Fig. 
[5l We find, in agreement with DMF T 32 i 33 and slave 
bosonSS calculations, that there is a critical value of hy- 
bridization strength (cf. Fig. [5j) below which the system 
becomes fully polarized (phase III) . For large value of hy- 
bridization, the strong ferromagnetic phase is destroyed 
because of the screening the f-electron magnetic moments 
partly by the c-electrons and partly by the itineracy of 
f-electrons. While a strong hybridization increases the 
magnetic interaction, it simultaneously weakens the mag- 
nitude of magnetic moments by the itineracy. When the 
hybridization parameter \V\ falls below the critical value 
system transforms in a discontinuous manner from the 
strong ferromagnetic phase III to the weak ferromagnetic 
phase. In the inset to Fig. [5] one can observe that the 
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FIG. 5: (Color online) The magnetization of c- and f- 
electrons, m c and m/, respectively, plotted as a function of 
hybridization parameter for selected numbers electron occu- 
pancy n e . The inset shows the corresponding dependence of 
the occupancy rif. 



occupancy rif simultaneously undergoes a sharp transi- 
tion as a function of the hybridization strength \V\. The 
limit rrif -» 1/2 and rif -> 1 (for \V\ < 0.15) can be 
thus called the localized moment limit, in which indeed 
V/e f « 1. 

The density-of-states (DOS) enhancement at the Fermi 
level is a direct measure of the effective mass enhance- 
ment of quasiparticles^— . As can be seen from Fig. [6j 
in the case of stable state II the density of states for the 
spin-down subband increases with the increasing the hy- 
bridization and that the spin- up part is absent, because 
the chemical potential is placed in hybridization gap for 
the spin-up subbands. 

The phenomenon such as the metamagnetic behav- 
ior of the magnetization curve, as well as the field- and 
spin-dependent effective masses have been observed in 
several heavy fermion compounds^"—. For selected pa- 
rameters, our model predicts that the metamagnetism 
appears in the crossover region in between the regimes 
II to III (see Fig. [7j). Namely, the principal feature of 
the heavy fermion state II (WFM phase) is its insensi- 
tivity to the strong applied magnetic field 40 . The only 
quantity which changes significantly is the chemical po- 
tential position. With increasing value of magnetic field, 
the chemical potential shifts and eventually leaves the 
quasiparticle-gap regime (cf. Fig. [8]), entering into the 
spin-up upper hybridized band and switching the system 
to phase III. This phase transition causes sudden increase 
of the magnetization and rearrangement of the quasipar- 
ticle subbands as can be seen in the insets to Fig. [51 
In the inset to Fig. one can see that define jump of 
rif in the applied magnetic field is observed while sys- 
tem evolves from weakly polarized (II) towards strongly 
polarized phase (III). The present result is consistent 
with the existing considerations, where additionally the 
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FIG. 6: (Color online) Density-of-states DOS enchancement 
at the Fermi level as a function of hybridization (a) and ap- 
plied field (b), illustrating the nature phase transition from 
phase II to III. The parameters are the same as those used 
above. The arrows indicate the spin orientations. 



Coulomb interaction between local and conducting elec- 
trons were taken into accoun t 41 ! 42 . The inclusion of the 
inter-orbital Coulomb interaction amplifies the sharpness 
of the transition. 

To understand better the behavior of strongly polar- 
ized HF system in phase III we analyze also the field 
dependence of the quasiparticle mass. As can be seen 
in the Fig. [9j the quasiparticle mass is suppressed with 
increasing magnetic field. It can be interpreted as a de- 
struction of the HF state. The magnetic field tends to or- 
der of both electrons and in this manner to suppress the 
antiferromagnetic Kondo coupling. The effective mass 
reduction at high fields is one of the characteristic fea- 
tures of HF— . In contrast to the paramagnetic ground 
state of PAM 45 , we observe a completely different be- 
havior of the effective-mass field dependence. The first 
difference is that spin-up mass is larger than the spin- 
down carrier mass due to the discussed above rearrange- 
ments of quasiparticle the subbands (see Fig. [8]), which 
is also consistent with the overall mean-field theory of 
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FIG. 7: (Color online) The magnetizations of m c and m/, 
plotted as a function of applied magnetic field. The metamag- 
netic (discontinuous) transition reflects the critical changeover 
when crosing the regime WFM— )>SFM. Inset: f-band occu- 
pancy change accompanying the metamagnetic transition. 

Wasserman et. al.— (who does not have the spin-split 
masses though). Secondly, both effective spin-dependent 
masses (a = ±1) decrease in the stronger applied field. 
The effective mass of quasiparticles depends also from 
the total number of electrons n e . In the case of the fixed 
hybridization strength the value of localized electrons n f 
is almost insensitive to the change in the total number 
of electrons, since in the region with an almost localized 
f-electrons (rif —> 1), changes in n/ are minor. Thus the 
n e dependence of effective quasiparticle mass in that re- 
gion is mainly due to change in number of conduction 
electrons n c . 



C. Application to Ce x L,aii- x ^G 

The obtained behavior can account for the field 
dependence of the quasiparticle effective mass of 
Ce^Lai-^B e 38143 in very strong fields. In this material, 
the dHvA signal is found to originate only from a single- 
spin-orientation Fermi surface sheet. However, it was dif- 
ficult to verify, from which spin channel it comes from. 
From the experimental observation^ 3 - it seems that spin 
up channel progressively dominates the dHvA signal as 
the concentration x is increased (i.e. the number n e in- 
creases). The experimental investigation claim that spin- 
up component may be dominant in this compound. In 
our approach we reproduce two important experimental 
features: (i) the field suppression of effective masses at 
high magnetic field, (ii) the non-monotonic x dependence 
of m* for the dominant spin component (cf. inset to Fig. 
[9j) and its small value for this heavy-fermion compound. 

One has to mention that Ce^Lai-^Be shows the Fermi 
T 2 - liquid behavior in low-temperature resistivity^ even 
though it exhibits a complicated phase diagram on the 




-1 -0.5 0.5 1 

CO— (J. 




-1 -0.5 0.5 1 

CO-JLl 



FIG. 8: (Color online) Spin resolved quasiparticle den- 
sity of states (DOS) in WFM (a) and SFM (b) phases. 
The upper (lower) plots corresponds respectively to the spin 
up(down) electrons. Insets: the band substructure for spe- 
cific parameters. The arrows label the spin orientation of the 
quasiparticles. The system is metallic in both WFM and SFM 
phases, but there is a change of spin polarization at Fermi en- 
ergy. 



magnetic field-temperature planed— . The discussed 
here effective mass suppression was observed in very 
strong fields, when the complicated quadrupolar phases 
discussed in^S have already disappeared. This is the rea- 
son why we can apply our result to this complicated sys- 
tem in strong fields. 



IV. SUMMARY AND OUTLOOK 

In the paper we have analyzed the effect of spin- 
split and applied-field dependent quasiparticle masses (or 
quasiparticle density of states) and their influence on se- 
lected magnetic properties. The results differ from those 
obtained for a single narrow band model. What is equally 
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FIG. 9: (Color online) A: A plot of the quasiparticle effective 
mass of the frequency versus concentration at 40 T (left 
of top figure A), the field dependence of the effective mass 
for different concentrations (right of top figure A). Figure A 
represents experimental results for CexLai-xB o 38 ^ 3 . Note 
that for A the parameters specified the mass enhancement 
is rather small. B: The field dependence of the quasiparticle 
effective masses for |V| = 0.3, remaining parameters are the 
same, as used before. The inset shows the quasiparticle ef- 
fective masses (in units of bare mass) and their dependence 
from the total number of electrons n e ; for applied magnetic 
field h = 0.03 and \V\ = 0.3. 



important, the heavy masses are strongly suppressed at 
the border of weakly- to strongly-ferromagnetic phase 
transition, as observed in Ce^Lai-^Be syste m 38 ! 43 . 

On formal side, we would like to note that we in- 
cluded large (but finite) U effects within the Ander- 
son lattice model, with both the hybridization and the 
Kondo interaction effects incorporated in a systematic 
manner. In result, the effective Hamiltonian can be 
termed Anderson-Kondo lattice model. The systematic 
aspect of our approach to solving the model relies on uti- 
lizing the Fukushima type of approach and generalizing 
it to a two-orbital situation, in which the constraints en- 
suring the statistical consistency of the results have been 
also included. In this manner, we can claim that we 
have provided a mutually consistent mean-field approach 
starting from the Gutzwiller-type projection for the cor- 
related f-electron states in this two-orbital model. The 



Kondo-lattice heavy-fermion fi uct uating-valence 
regime: (Anderson-Kondo) regime: 

Hamiltonian (5) regime: Hamiltonian (4) 
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i 1 1 i i i i i i 
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FIG. 10: (Color online) Schematic representation of various 
regimes as denned by effective model Hamiltonians. The na- 
ture of the borders between the regions is not sharply defined. 
This classification is particularly relevant for various versions 
of the so-called renormalized mean- field analysis. 



Appendices [A] to O provide some of the technical details 
of our approach. 

Generally, it comes out from our results (as well as from 
the earlier results^£) that for \V\ < \V min \ (~ 0.25 for 
ef = —0.75) the f-electrons are almost localized (rtf — )• 1 
c.f. inset, Fig. [5j see also Fig. 13(c)). Well below the criti- 
cal value |Knin| we can say that the Kondo-lattice Hamil- 
tonian (5) is applicable. In the heavy-fermion regime, we 



have \ Vrr, 



< 



V\ < \ef\ and the Anderson-Kondo lattice 
model in the forms (j3]) or (4) should be used, whereas 
V/ef > 1 we enter a fluctuating valence regime, where 
the Falicov-Kimball term may also be important to cap- 
ture the large change in n/ : n c . The discussed sub- 
division into regimes is illustrated schematically in Fig. 
[TOl Obviously, in all this consideration we assume that 
|V|/(C/ + e/ )«l. 

Finally, the approach is an effective (almost localized) 
Fermi (non-Landau) picture with renormalized charac- 
teristics nonlinear effective fields (A/, A/ m , A cm ), spin- 
dependent effective masses and adjustable chemical po- 
tential in each phase. All the parameters are determined 
in a consistent manner by determining first the effective 
Gibbs energy functional which plays role of Landau func- 
tional with number of the mean fields (order parameters). 
It would be natural to apply the scheme to the analysis of 
the paired (superconducting) states. This is planed as the 
next step. Also, a coexistence of (antiferro) magnetism 
and superconductivity should be analyzed, since it is ob- 
served in quite few heavy-fermion compounds^—. 
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Appendix A: Role of the interorbital 
(Falicov-Kimball) term 

The Falicov-Kimball term taken in the intraatomic 
form 



(Al) 



enters in two ways in the starting effective Hamiltonian 
([3]). First, the denominator of the third term should be 
replaced by U — U c f + ey, so that the effective Kondo 
coupling constant increases, since it has then the form 
Ji m = 2\Vi m \ 2 /(U + U c f + ef). Second, after the canon- 
ical transformation of (pQ) leading to ([3]) there remains 
the projected part of repulsive Coulomb interaction term 
Ucf^^myVifhm which, when treated in the simplest 
Hartree-Fock manner, reduces to 



cf 2^ V if n ™ 

<im> 



U c f (rifh m + n c % - ricnf). 



J cf A 

<im> 



(A2) 

In this manner, this term first of all shifts the energy / 
and c levels. Since such shifts are already included in (|9j) 
via the terms ~ A/ and A, respectively, to be calculated 
self-consistently anyway, we can disregard them. This is 
particularly so in the interesting us heavy- fermion- liquid 
regime (Oif) = 1 (jif — >> 1), and then 



<im> 



<im> 



[pif - i)h m + u C f 



<im> 



(A3) 

In effect, the first term is relatively small and hence can 
be decoupled in the Hartee-Fock manner. Such reasoning 
justifies qualitatively the fact of neglecting the Falicov- 
Kimball term in our mean-field-like approach, at least in 
the heavy- fermion limit. 



where projector V can be presented in the following form 

-Pi = Aj^^l-iV^t) 

= (1 - N it )(l - Nil) + \A^(1 " N it ) 
+V / V^t( 1 -^i)- (B3) 
The uncorrelated average of the norm is 



(Pt)o 



2\ _ (1 -^i/tX 1 ~ n ifi) 



f 



(B4) 



Using the Wick theorem we can easily calculate the fol- 
lowing averages as 

(ViN if& (l - N ia )c ma c ma Ti) Q 
= (\i<jNi<j(l - A r i cr )c] no .c ma )o 
= A^ [riifa(l ~ nif a )n iC a + ^/aTwTw] , (B5) 



Vt C ma ficr fi a C"m<jPi ) 



( V ' ^ia^ia c lna ficr fi a c ma ) 

VaA^wtL- ( B6 ) 



Finally, 



{Ki) = -UifaUica ~ — 



Mfa 



l ifcr 



r Yimcr r Yi; 



yJ{l-n ifo -)(l-n ifa ) 



(B7) 



Appendix C: Density of states and effective- mass of 
renormalized quasiparticles 

Using quasiparticle energy spectrum (JTTJ we can eval- 
uate the density of states by performing the sum 



Appendix B: The Fukushima local-constraint 
scheme of evaluating the expectation values 

In order to provide an insight into the procedure of cal- 
culating the expectation values in the Fukushima local- 
constraint scheme, we present here the detailed eval- 
uations of the average representing antiferromagnetic 



Kondo coupling K = Si ■ s n 



from Eq. 



4 1 4 

(4). After performing simple algebraic manipulations we 
can present this term in the following form: 



Ki 



-N if& (l - Nif^C^Cma ~ C m „fiaficm*- (Bl) 

According to (j6j) we can present expectation value of 
above operator as 



(Ki) 



(VN i9 {l-N i(T )cl na c m(J 'P) Q 
(V 2 )o 

(Pclnafi*fLcma'P)o 

(P 2 )o 



(B2) 



(CI) 



This can be carried out numerically by representing delta 
function as the Gaussian function with a small finite 
width. We may also achieve an analytical form for DOS, 
as follows. The matrix of the coefficients (jC2|) from the 
equation ([TO]) can be shifted^ so that obtained matrix 



has a triangular form 

JCx + (\ + \ f + \ fm a + ii)I 



6c - b -T 
-T a 



(C2) 



(C3) 



where b — —A/ — cr(A/ m — A cm ) and £k = —2r](cosk x + 
cos k y ). The eigenvalues of the above matrix are as fol- 
lows 



& - b + n V(& - b) 2 + 4rj) . 



(C4) 
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Then, it easy to see that the sum 

D e ( W ) = ^^-eL) 



kern 



(C5) 



corresponds to the true DOS ()C1|) shifted by the value 
A + A/ + A/ m a + fi To simplify the summation (jC5|) we 
use properties of delta function so that D e (uo) can be 
rewritten in the form 



kcr 



h y + ^5{^-u:{^-b)-4). 

(C6) 

In the thermodynamic limit the energy spectrum is ex- 
tremely dense so that it can be presented as smooth func- 
tion and summation can be substituted by integral so 

2_ 2 

that — >> cj. Using this we can express (£k— b) —> " J° 
so that 



D e (uj) 



kcr x ' x 



1-^ 



(C7) 

For the noninteracting electrons density of states can 
be written as the product of a line-shape function and 
a heaviside function, which ensures that the density of 
states vanishes outside the band: D(uo) = g(uj)S(W 2 — 
4cj 2 ), where W is the width of a band. In the similar way 
manner Eq. (|C7|) can be expressed in the form: 



f(uj) =0 and are as follows 



^3cr 
^4<T 



V(6/2 + W/2) 2 +r2 



■y/{b/2-W/2Y 



y/(b/2 + W/2) 2 



V# - ^/2) 2 



- W74 - 6/ 2 » 
2 -L- W/A-b/2, 
*-W/4-b/2, 
J+W/A-b/2. (C9) 



The DOS has non-zero values only in two intervals 
(o;i cr ,cj2cr) and (c^ , ) . The resulting expression for 
D e (uj) is finally as follows: 



D *M = ^E( 1 + ^) e ^- ' ( C1 °) 



We define the effective mass enhancement of quasipar- 
ticles in the standard manner 



m 



dE\ 



kcr 



de 



(Cll) 



Using the same substitution as above (i.e., — >• gj), 
we can express the starting energy of free electrons by 

2 2 f 

so that finally the mass enhancement 



at the Fermi energy is 



(C8) 



d f uo 



uj=0 



where f{u) = b + (^1 - ^J. The edges of DOS for xhis expression is used for numerical calculations in Sec. 
interacting DOS are given by the solution of equation HIB 
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